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Abstract 

Let X = {X{p),p ^ M} be a centered Gaussian random field, where M is a smooth 
Riemannian manifold. For a suitable compact subset D G M, we obtain the approx¬ 
imations to excursion probability Pjsuppg^, X(p) > u}, as rt —>■ oo, for two cases: (i) 

X is smooth and isotropic; (ii) X is non-smooth and locally isotropic. For case (i), the 
expected Euler characteristic approximation is formulated explicitly; while for case (ii), it 
is shown that the asymptotics is similar to Pickands’ approximation on Euclidean space 
which involves Pickands’ constant and the volume of D. These extend the results in 
from sphere to general Riemannian manifolds. 
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1 Introduction 


For a Gaussian random field X = {X{p),p £ M} living on some parameter space M, there 
is rich literature on studying the excursion probability Pjsuppg^, X(p) > u} as u ^ oo, 
where D is a suitable compact subset of M [l|, 0,0, [l^ . Traditional research usually focuses 
on Gaussian random fields living on Euclidean space. However, due to recent development 
in statistics and related fields such as geoscience, astronomy and neuroscience, Gaussian 
random fields parameterized on manifolds are becoming more and more useful for modelling 
many research objects. Motivated by both real applications and theoretical development 
0,0, BBS, 0], it is valuable to extend the study of excursion probability from Euclidean 
space to manifolds. 


When X is smooth and M is a Riemannian manifold, Taylor et al. 17j (see also Theorem 
14.3.3 in B) showed that the excursion probability can be approximated by the expected 
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Euler characteristic of the excursion set of X exceeding u. This verifies the so called “Euler 
Characteristic Heuristic”. The expected Euler characteristic can be formulated by a nice but 
implicit form (see [l^ or Theorems 12.4.1 and 12.4.2 in j^), which involves Lipschitz-Killing 
curvatures of the excursion set. Since the Lipschitz-Killing curvatures depend on both the 
field X and the parameter space M, it is usually hard to obtain their explicit expressions. 
In this paper, we consider X to be isotropic on M in the sense that its covariance function 
can be written as C{p, q) = p{d‘j^{p, q)) for any p,q ^ M, where p is a real function and cIm 
is the geodesic distance on M. It is then showed that the Lipschitz-Killing curvatures of the 
excursion set can be simplified as depending only on p^(0) and the geometry of D, see Theorem 
O below. This leads to a relatively simple formula for the expected Euler characteristic of 
the excursion set, which approximates the excursion probability with a super-exponentially 
small error as n —>■ oo (see Theorem 12.31) . 

We also study another case when X is non-smooth and M is a Riemannian manifold. In 
particular, we assume the covariance function of X satisfies the locally isotropic condition 
(13.11) below. Notice that since X is non-smooth, the expected Euler characteristic approxi¬ 
mation is no longer applicable. Instead, we shall employ arguments by charts to study the 


corresponding field under Euclidean local coordinates and then apply the result in to find 
an approximation to the excursion probability. As a related reference, we have found that 
Mikhaleva and Piterbarg 11[ studied the excursion probability over a manifold subset for 
Gaussian random fields parameterized on Euclidean space. Their theorems can be applied 
to obtain results similar to Theorem 13.41 in this paper if X is the restriction on M of some 
Gaussian field defined on Euclidean space of higher dimension. However, it is not clear if 
every locally isotropic Gaussian field on M can be regarded as a restriction of some Gaussian 
field on Euclidean space satisfying the local property in (4) in 11[. To see this, we refer to 


explanations in for the case when M is an A^-dimensional sphere (the case when M is an 
arbitrary Riemannian manifold will make it even more complicated). On the other hand, in 
many real applications, it is natural to model a Gaussian field on M directly and then explore 
its local properties on M but not its extension on Euclidean space of higher dimension. In 
such sense, condition (13.Ih is more natural and useful compared with (4) in lll |. 


The results in this paper extend those in 


from sphere to general Riemannian man- 

relies on 


ifolds. The study for excursion probabilities of Gaussian fields on sphere in 
specific spectral representation of the covariance [l^ and the spherical coordinates. To deal 
with Gaussian fields on general Riemannian manifolds, we explore directly the connection 
between the Riemannian metric and the (locally) isotropic property of the field to obtain the 
asymptotics of excursion probabilities (see Theorem 13.41 and Gorollarv 13.51 below). This new 
approach is more general and powerful. 
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2 Smooth Isotropic Gaussian Random Fields on Manifolds 


2.1 Preliminaries on Riemannian Manifolds 

We first introduce some notations and definitions on Riemannian manifolds. Let (M, g) be 
an A^-dimensional smooth Riemannian manifold, where g is the Riemannian metric. Then 
for each p G M, there is an inner product 


gp : TpM x TpM —>■ M 

{^p: ^p) 9p(Api ^p)i 

where TpM is the tangent space of M at p. To simplify the notation, we denote by {^p, ap) 
the inner product gp{^p,(Tp). In particular, for a coordinate system with bisis {-^\p}i<i,j<N, 
we let 


(A. 

d 


/A 

d 

[dxi 

p’ dxj 

J- 

\dxi 

p’ dxj 


which are called the components of s' at p under this coordinate system. Denote the N x 
N matrix by G{p) = {gij{p))i<i,j<Ni which is symmetric and positive dehnite due to the 
definition of g. For a real-valued smooth function / on M, let fi = and /jj = 

Notice that, the Riemannian metric g itself is not a true metric on M. However, it does 
induce a metric on M in the following way. Firstly, we can define the length of a differentiable 
curve 7 : [a, 6 ] —>■ M by 


L(7) = / 


Then the distance on M induced by the Riemannian metric g, denoted by dMi', •)) is dehned 


by 


dM{p,q) =infL{'y), p,qeM, 


where the inhmum extends over all differentiable curves 7 beginning at p & M and ending at 
qe M. 


2.2 Isotropic Gaussian Random Fields on Riemannian Manifolds 

Let X = {X{p),p G M} be a smooth centered Gaussian random field parameterized over 
an Wdimensional Riemannian manifold {M,g). We define X to be isotropic over M if the 
covariance function satisfies 


C{p,q) = p{d\j{p,q)) , p,qeM, ( 2 . 1 ) 

where p is a real function on [0,oo). Relation (12.ip implies that C{p,q) only depends on 
dM{p,q) and hence behaves isotropically over M. Besides the well-known Euclidean case, we 
refer to [l^, and respectively for isotropic Gaussian fields on sphere and hyperbolic space 
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satisfying m- For general Riemannian manifolds, isotropic G^ssian fields and covariances 
can be constructed via spectral representations; see 2{j, [ijl, ll 


As discussed in (see Corollaries 11.3.3 and 11.3.5 therein), in order to apply the 
Kac-Rice metatheorem and the Morse theory to establish the expected Euler characteristic 
approximation, we need the following smoothness and regularity conditions for X. 

(HI). A(-) G C'^{M) almost surely and there exist positive constants K and rj such that 

E{X,j{p)-Xij{q)f <K\log{dM{p,q))\-^^^^\ yp,qeM, i,j = l,...,N. 


(H2). For every p G M, the Gaussian vector {Xi{p), Xjk{p),l < i < N,1 < j < k < N) is 
non-degenerate. 


In addition, to make the compact subset D <Z M suitable, we shall assume that H is a 
regular stratified, locally convex submanifold. We refer to pages 198 and 189 in for 
rigorous dehnitions of “regular stratified” and “locally convex”. Roughly speaking, we may 
let D be piecewise smooth and convex. 

Next, we introduce the Lipschitz-Killing curvatures (or “intrinsic volumes”) of D. For 
j = 0,, dim(Zl), the j-th Lipschitz-Killing curvature Cj{D) can be viewed as the measure 
of the j-dimensional boundary of D. Generally, if Zl is a Z-dimensional, regular stratified, 
locally convex manifold, then Cq{D) is the Euler characteristic, Ck{D) is the volume and 
Ck-i{D) is half of the surface area. For the other Cj{T), 1 < j < /c — 2, we can apply 
Steiner’s formula (see page 142 in Q]) to find their values. We refer to Q] for the general 
definition of Cj{D) and specific examples on evaluating them explicitly as well. 

Let ^Au(X, D)) be the Euler characteristic of excursion set Au{X, D) = {x ^ D : X{p) > 
u} (cf. [^). Denote by Hj{x) the Hermite polynomial of degree j, that is, 

«,(x) = 

We obtain the following formula for x(A„(X, D)) when X is smooth and isotropic. 


Theorem 2.1 Let X = {X{p),p G M} 6e o centered, unit-variance, isotropic Gaussian 
random field satisfying (El]), (HI) and (H2), where (M, g) is a smooth Riemannian manifold. 
Then, for a compact, regular stratified, locally convex submanifold D C M , 

dim(D) 

E{x(vl.(X,D))}= (-2p'(0))''/'£,(D)/3,(u), (2.2) 

i=o 


where Cj{D) are the Lipschitz-Killing curvatures of D and 


/3j{u) = 


(27r) e 1‘^dx 

(27r)-(-?+i)/2iZj_i(n)e-“"/2 


*/j = 0, 

if 2 = l,...,dim(D). 


4 



Proof By Theorem 12.4.2 in [^, 

dim(D) 

¥.{x{Au{X,D))}= Y. Cf{D)P,{u), 
j=o 

where C^{D) are Lipschitz-Killing curvatures of D under the Riemannian metric induced 
by X. Specifically (see page 305 in [^), 

dpoi^poi ^po) ■ {^poX)'\ = ipoO'poC{p^Q)\p=q=p(n Vpo £ XI, 

where Cpo)<7po ^ Tp^M, the tangent spaces of M at po- Therefore, to prove ()2.2p . we only 
need to show Cf{D) = {—2p'{^)y^‘^Cj{D) for j = 0,... ,dim(T)), where Cj{D) are original 
Lipschitz-Killing curvatures of D under the Riemannian metric g. 

We may choose two geodesic curves on M, say 7 : [0, <5] —>■ M and r : [0, (5] —M, such that 
7 (t) = exppp(t^p(,) and r(s) = expp|^(s(Tpp), where exp is the usual exponential mapping on 
M. Here 5 is small enough so that the squared distance function d\p{'^{t),T{s)), t,s G [0,(5], 
is smooth. 

Let t G (0, 5) be fixed. For each s G [0, (5], consider the minimizing geodesic ps ■ [0,1] —M 
joining 'y(t) to r(s), that is, 7s(0) = 7(1) and 7s(l) = t(s). Dehne the variation f{s,x) = 
7s (x), where x G [0,1] and s G [0,(5], and dehne the energy function 

lo ^ dlp{-/{t),T{s)), 

where the last equality is due to the fact that 7s(') is geodesic for each hxed s. By the hrst 
variation formula for energy (see Proposition 2.4 in Chapter 9 of 0]), together with again the 
fact that 7s(-) is geodesic, we obtain 

^c^M(7(i),T(s)) = E'{s) = 2(T'(s),exp"(^^)(r(s))) = 2((Tpo, exp“(^^^ (r(s))), 

where exp“^^^(r(s)) denotes the vector ( G such that t { s ) = exp.^^^)(C). In particular, 

when s = 0, exp“^^^(T(s)) = —t^pQ- Therefore, 

'^d,]yf{'y(t),T{s))\s=o = — 2t((Tp(j, ^pp). (2-3) 

It follows from (12.3p that 

d d 

^PO^PoE{Pi(l)\p=q=po — ('^m( 7(^))'^(■s))) |i=s=0 

d 

= -2^ {dli{l{t),po))] [ 4=0 

= ~2p (0)((TpQ, ^po)- 
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Hence the Riemannian metric induced by X is given by 

dpoi^poi ^po) = (0)(i^pQ, ct^q), \/po G M, ^poj'^po ^ Tp^M. 

By the definition of Lipschitz-Killing curvatures, one has Cj{D) = {—2p'{Q)y^‘^ Cj{D), yield¬ 
ing the desired result. □ 


Remark 2.2 The following are some remarks. 

• Let M = that is, {X{p),p G M'^} is an isotropic Gaussian random field on Eu¬ 
clidean space. Then, since Var(Xj(p)) = —2p'{0), we see that (12.2.19) in gives the 
same result as Theorem EH 


Let M be the A^-dimensional unit sphere S'^. For simplicity (the general case of co- 


variance 


[^, 15] can be handled similarly), we assume the covariance of the isotropic 


Gaussian field on sphere has the form 

OO 

C{p,q) = '^bn{p,q)'^N+i, p,qeS^, 


n=0 


where > 0 and (•,-)]gjv+i denotes the usual inner product in Notice that 

ds^iP^q) represents the spherical distance (or the greatest circle distance), therefore. 


OO 

C{p, q) = '^bn (cos dgiv (x, y)T = p {dl^ {x, y)) , 

n=0 


where p{r) = ^ Then —2p'{0) = which is the same as 

the parameter C m. [5|. This indicates that the result for isotropic Gaussian fields on 
sphere in (especially Lemma 3.5 therein) is a special case of Theorem 12.11 


• It can be seen from the proof of Theorem 12.11 that, under an orthonormal basis of the 
Riemannian manifold {M,g), say {-^\p}i<i,j<N■, one has 

¥.{X,{p)Xj{p)} = -2p'(0)5,,-, Vp G M, 1 < ij < N, 


where 6ij is the Kronecker delta function. 

Applying Theorem l2.1l and Theorem 14.3.3 in Adler and Taylor j^, we obtain immediately 
the following approximation for the excursion probability. 

Theorem 2.3 Let the conditions in Theorem \2. 1\ hold. Then, following the notation therein, 
there exists oq > 0 such that as u ^ oo, 

( 1 dim(D) 

P supA(p)>u = (-2p'(0))^/'£,(Il)/3,(u) + o(e-“°"/2^. 

J j=o 
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The approximation provided in Theorem 12.31 is very accurate since the error is expo¬ 
nentially smaller than the expected Euler characteristic. In contrast, the approximation for 
non-smooth Gaussian fields only gives the leading term (i.e., the largest order of polynomials 
of u) of the expected Euler characteristic, which can be seen in 13|] and the results in the 
section below. Therefore, in many statistical applications, especially when using excursion 
probability for computing p-values, the expected Euler characteristic approximation is much 
more popular. 


3 Locally Isotropic Gaussian Random Fields on Manifolds 

As before, let (M, g) be a smooth Riemannian manifold of dimension N. In this section, 
we consider a centered Gaussian random field X = {X{p),p G M} such that the covariance 
function satisfies 


C{p,q) = 1 - cdM{p,q){l + o{l)) as dM{p,q) ^ 0, 


(3.1) 


where c > 0 and a G (0,2]. If X is smooth, such as satisfying condition (HI), then a = 2. 
For a G (0,2), the field X can be checked to be non-smooth. Covariances satisfying (|3.ip 
behave isotropically in a local sense, therefore we call them (or the Gaussian fields) locally 
isotropic. In particular, if the covariance is isotropic with the representation ()2.ip . then ()3.ip 
usually holds. For the Euclidean case, i.e., M = there is rich literature on stud ying the 


12 
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3B). 


asymptotics of excursion probability of Gaussian fields satisfying ()3.ip (cf. 

Recently, Cheng and Xiao B studied Gaussian fields on sphere satisfying (13.11) by applying 
the spherical coordinates. Here, we shall use the Riemannian metric itself to show that the 
locally isotropic property can be generalized as (|3.1I) for studying the asymptotics of excursion 
probability of Gaussian fields on Riemannian manifolds. 

The main approach is using the argument by (coordinate) charts to transform condition 
(|3.1h to an equivalent condition in Euclidean space (see Lemma 13.11 below). We first recall 
some elementary terms for manifolds, which can be found in many textbooks such as B, Bj- 
A chart on M is a pair {U,(p), where U C M is open and p : U ^ C is a 

homeomorphism. If a collection of charts, say gives a covering of M, i.e., Uig/I/j = 

M, then it is call an atlas for M. The component functions xi,... ,xn of a chart {U,(p), 
defined by (p{p) = {xi{p),... ,xm{p)), are called the local coordinates on U. 

For two functions / and h, we say f{t) h{t) ast^ to £ [- 00 , -l-oo] if limt^t^ f{t)/g{t) = 
1. Denote by || • || the usual Euclidean norm. 

Lemma 3.1 Letp, q £ U, where p is fixed and {U, p) is a chart of M. Then as dM{p, q) —^ 0, 

N 

dhip^q) ~ 9ij{p)[xi{q) - Xi{p)][xj{q) - Xj{p)] = \\G^^^{p) • ip{q) - t{p))\\^ , 
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where G^/‘^{p) is the square root of the positive definite matrix G{p) = igij{p))i<i,j<N defined 
in Seetion [23 


Proof Let 7 : [0, 1] —)• M be a smooth curve on M such that 7 ( 0 ) = p and 7 ([ 0 , 1]) C U. 
For J > 0, denote by L( 7 ([ 0 ,( 5 ])) the length of the segment between 7 ( 0 ) and 7 (( 5 ). Since 


N 


i'{t) = Y. 


2=1 


dxijjjt)) d 
dt dxi ’ 


we obtain that, as <5 —)• 0, 


^( 7 ([ 0 ,< 5 ])) = ^ ^g^(t){i{t),i{t))dt 


1/2 




dt 


dt 


N 


1/2 


This implies 


N 


Y 9ijip)[xihiS)) - Xi{p)][xj{'y{5)) - Xj{p)] 


1/2 


dM{p,q) ~ I Y 9ijip)[xi{q) - Xi{p)][xj{q) - Xj{p)] | = \\G^^^{p) ■ (ipiq) - p{p))\\- 

77=1 


□ 


Remark 3.2 Let M be the A^-dimensional unit sphere S^. Following the notation of 
spherical coordinates in [^, we have 


N-l 


G{p) = diag 1, (sin 61 )“^, • • •, (sin 6 i) 


2=1 


Therefore, Lemma 2.1 in is a special case of Lemma l3.11 
Denote by Ha,N the usual Pickands’ constant in R^, that is 


Ha,N = lim K 


-N 


K^oo 




sup Z{s) >u \ du, 
6G[0,K]^ I 


where {Z(s) : s G [0,oo)^} is a Gaussian random field such that 
E(Z(s)) = -||s||“, Coy{Z{s), Z{v)) = ||s||“ + 


V — s — V 





Let T C be a bounded A^-dimensional Jordan measurable set (i.e., the boundary of 
T has A^-dimensional Lebesgue measure 0). Let Y = {Y{t),t E be a centered Gaussian 
field with covariance function Cy such that for every fixed t ^ T, 


Cyit, s) = 1 - ||iJ(t)(s - t)||“(l + o(l)) as ||s - t\\ 0, 


(3.2) 


where B{t) is a non-degenerate N x N matrix and a E (0, 2]. 

We will make use of the following result, which is_a direct consequence of Theorem 2.1 in 
(see also Theorem 2.2 in [^) and Lemma 2.3 in 


Theorem 3.3 Let T C be a hounded N-dimensional Jordan measurable set. Suppose the 
centered Gaussian random field {Y{t),t E satisfies condition h3.S\) . Then as u ^ oo, 

supy(t) > |det(i3(t))|dy -|-o(l)). 


Here and in the sequel, ^'(m) = (27r) e '"'^1‘^dv. 


Now we are ready to prove our main result of this section. 

Theorem 3.4 Let X = {X(p),p E M} be a centered Gaussian random field satisfying con¬ 
dition (M), where {M,g) is an N-dimensional smooth Riemannian manifold. Let D <Z M 
he an N-dimensional smooth compact submanifold on M. Then as u ^ oo, 


sup X{p) > u 
p&D 


+ o(l)), 


where Vol(iJ) denotes the volume of D and c > 0 is the constant in 113.1\) . 


Proof Let {Ui,ipi)i^i be an atlas of M. Since D is compact, it has a finite covering 
{Ui, (pi)i^iQ, where Iq C I is a finite index set. For i E Iq, define X : ^pfiUi) C by 

X = X o ip~^. By Lemma [3Tl the covariance of X, denoted by : (pi{Ui) x ipi{Ui) C 

X — )• M, satisfies the following property: for every fixed (pi{p) E (pi{Ui), as \\g:>i{q) — 

<^i(p)ll -S' 0, 

C{Ti{p),Ti{Q)) = 1 -c||(G^/2 o(^-i)((^i(p)) • {cpfiq) - 99i(p))||"(l -h o(l)) 

= 1 - o ip-^){ipi{p)) ■ {iffiq) - ipi{p))\fi{l + o(l)). 

We make a decomposition D = such that any neighboring Dj and Dji belong to a 
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same chart. Suppose Dj C Ui, applying Theorem 13.31 yields 


sup X(p) > u > = P < sup X{x) > u 

peDj I [xe<pi(Dj) 




'tPiiDj) 


det o (p^ (x)^ dxj //a,iV'W^^^“'I'(n)(l + o(l)) 

|det ((Go (^r^) (x)) ) c^/"Fa,Aru^^/“T(u)(l + o(l)) 


(3.3) 


= Vo1(Z?,)c'^/“F„,^u2^/“T(u)( 1 + o(l)), 

where the last line is due to the definition of volume on (M,g). 
It is obvious that 



On the other hand, by the Bonferroni inequality, 


(3.4) 


p < sup x{p) > u> > ] 

J i=i 


sup X{p) > u 

. per>i 


E 


P < sup X{p) > u, sup X{q) > u 

I p&Dj qGD,, 


(3.5) 


If Dj and Dj/ are not adjacent, i.e., dM{Dj,Dji) > 0, then it follows from the Borell-TIS 
inequality (cf. Theorem 2.1.1 in [^) that 

P { suppgo. X{p)>u, supgg£, X{q)>u\ 

lim —- — 

u—>-oo (^u) 

IPlsuPpe/j gen.,[^(P) + ^(9)]/2 > 

< lim —^--- - = 0. 

u^oo 

If Dj and Dj/ are adjacent, then by our assumption, there exists a chart (JJi,pi) containing 
both Dj and Dj/. Since 


P < sup X{p) > u, sup X{q) > u 

I PGDj q€D./ 


sup X(p) > u > + P < sup X{q) >u> — 


sup X{p) >u 


J \^q£Dj/ J [^pGDjUDj/ J 

applying (13.3p to the last three terms respectively yields that the joint excursion probability 
above is o(tt^'^/“'I'(tt)). Therefore, the last sum in (|3.5p is o(u^'^/“'I'(u)). 
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Now, combining (13.311 with (I3.4p and (13.5p . we obtain 

p|supX(p) >u \ = f" Vol{Dj)c^/'^ Ho, + o{l)) 

U 

= Vol(L>)c^/“/7«,7vu2^/"^(n)(l + o(l)). 

□ 

In Theorem 13.41 the submanifold D <Z M \s assumed to be of the same dimension as M. 
This is because under such assumption we can apply the Riemannian metric g on M directly. 
If dim(Z)) < N, we may turn to the restriction of g on the tangent space of D, which can 
be regarded as the Riemannian metric on D induced from {M,g). We denote such induced 
Riemannian metric on D by g. Then we have following result as an extension of Theorem 


Corollary 3.5 Let X = {X{p),p G M} be a centered Gaussian random field satisfying 
condition USA]) , where {M,g) is an N-dimensional smooth Riemannian manifold. Let D G M 
be a k-dimensional smooth compact submanifold on M, where 0 < A; < A^. Then as u ^ oo, 


P 


sup X{p) > u 
p£D 




where Xol{D) denotes the volume of D under the induced Riemannian metric g and c > 0 is 
the constant in (Ejp. 


Proof Let (Vi, be a finite atlas of D. For i G Jq, defnie X : (fiiVi) C —>• 

by X = X o Similarly to Lemma l3.ll we have that the covariance of X, denoted by 
: (j)i{Vi) X (j)i{Vi) C X ^ M, satishes the following property: for every fixed 
4>i{p) e 0i(Vi), as \\4>i{q) - Mp)\\ 

C{Mp),Mq)) = 1 - o (p-^){(pi{p)) ■ {(pi{q) - (pi{p))P{l + o{l)) 

= 1 - ||cV“(gV2 o cf-^){Mp)) • (Mq) - Mp))ra + 0(1)), 

where G{p) = {gij{p))i<ij<k- The desired result then follows from similar arguments in the 
proof of Theorem 13.41 □ 

Remark 3.6 If X = {X{p),p G M} is isotropic such that conditions (HI) and (H2) are 
fulfilled, then the covariance C{p,q) = p [d\j{p,q)) satisfies condition (13. ip with a = 2 and 
c = —p'{Q). Since Pickands’ constant H 2 ^n = one can check that the approximation 

in Theorem 13.41 only provides the leading term of the approximation in Theorem 12.31 This 
is because. Theorem 13.41 does not take into account the boundary effect of D. Therefore, 
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the expected Euler characteristic approximation in Theorem 12.31 is much more accurate for 
smooth isotropic Gaussian fields. However, Theorem 13.41 is still very useful since it provides 
an approximation to non-smooth locally isotropic Gaussian fields, for which Theorem 12.31 is 
not applicable due to lack of smoothness. 

In this paper, for the non-smooth case, we only study covariances with a simple local 
structure m- However, similarly to the Euclidean case (l^ . one can also investigate the 
asymptotics for Gaussian covariances with more general local structures, such as letting the 
constants c and a in (j3.1l) be functions c{p) and a{p) depending on p Q]. We leave these 
problems for future research. 
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